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We investigate a C/(l) lattice chiral gauge theory (LxGT) with domain wall fermions and compact 
gauge fixing. In the reduced model limit, our perturbative and numerical investigations show that 
there exist no extra mirror chiral modes. The longitudinal gauge degrees of freedom have no effect 
on the free domain wall fermion spectrum consisting of opposite chiral modes at the domain wall 
and at the anti-domain wall which have an exponentially damped overlap. 
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I. INTRODUCTION 

Lattice regularization of chiral gauge theories has remained a long standing problem of nonperturbative investigation 
of quantum field theory. Lack of chiral gauge invariance in LxGT proposals is responsible for the longitudinal gauge 
I degrees of freedom (dof) coupling to fermionic dof and eventually spoiling the chiral nature of the theory. The well- 
known example is the Smit-Swift proposal of LxGT 0]. Although in a recent development using a Dirac operator 
that satisfies the Ginsparg- Wilson relation, it was possible to formulate a LxGT without violating gauge-invariance 
or locality ||^, an explicit model for nonperturbative numerical studies is still not available. 

In this paper, we follow the gauge fixing approach to LxGT |^ . The obvious remedy to control the longitudinal gauge 
dof is to gauge fix with a target theory in mind. The Roma proposal Q| involving gauge fixing passed perturbative 

• tests but does not address the problem of gauge fixing of compact gauge fields and the associated problem of lattice 
^ I artifact Gribov copies. The formal problem is that for compact gauge fixing a BRST- invariant partition function 

0^ I as well as (unnormalized) expectation values of BRST invariant operators vanish as a consequence of lattice Gribov 
copies Shamir and Golterman |^ have proposed to keep the gauge fixing part of the action BRST noninvariant 
and tunc countcrterms to recover BRST in the continuum. In their formalism, the continuum limit is to be taken 
from within the broken ferromagnetic (FM) phase approaching another broken phase which is called ferromagnetic 
I directional (FMD) phase, with the mass of the gauge field vanishing at the FM-FMD transition. This was tried 

• out in a U{1) Smit-Swift model and so far the results show that in the pure gauge sector, QED is recovered in the 
continuum limit and in the reduced model limit (to be defined below) free chiral fermions in the appropriate chiral 

d . representation are obtained ||^. Tuning with countcrterms has also not posed any practical problem, actually very 
I ' little tuning is necessary. Efforts are currently underway to extend this gauge fixing proposal to include nonabclian 
gauge groups g. 

~5 ' Without gauge fixing the longitudinal gauge dof, which are radially frozen scalar fields, are rough and nonpertur- 
' bative even if the transverse gauge coupling may be weak (this is because with the standard lattice measure, each 
point on the gauge orbit has equal weight). The theory in the continuum limit, taken at the transition between 
the broken symmetry ferromagnetic (FM) phase and the symmetric paramagnetic (PM) phase, displays undesired 
nonperturbative effects of the scalar-fermion coupling that usually spells disaster for the chiral theory. The job of the 
' gauge fixing is to introduce a new continuous phase transition, from the FM phase to a new broken symmetry phase 
(FMD), at which the gauge symmetry is recovered and at the same time the gauge fields become smooth. 

The problem can be cleanly studied in the reduced model as explained in the following. When one gauge transforms 
a gauge non-invariant theory, one picks up the longitudinal gauge degrees of freedom (radially frozen scalars) explicitly 
in the action. The reduced model is then obtained by making the lattice gauge field unity for all links, i.e., by switching 
off the transverse gauge coupling. The action becomes that of a chiral Yukawa theory with interaction between the 
fermions and the longitudinal gauge dof. The reduced model would have a phase structure similar to the full theory, 
e.g., the gauge fixed theory in the reduced limit will have a FM-FMD transition in addition to the FM-PM transition. 
Now for the gauge fixing proposal to work, the scalars need to decouple from the fermions at the FM-FMD transition 
leaving the fermions free in the appropriate chiral representation. Passing the reduced model test is an important 
first step for any LxGT proposal that breaks gauge invariance. 

In the reduced model derived from the gauge fixed theory the scalar fields become smooth and expandable in a 
perturbative series as 1 -I- ©(coupling constant) at the FM-FMD transition. If continuum limit can be taken near 
the point in the coupling parameter space around which this perturbative expansion is defined, the scalar fields will 



X 



*Permanent Address: Dept. of Physics, NND College, Calcutta 700 092 



1 



decouple from the theory. The parameterization of the gauge fixing action turns out to be a good one, because this 
continuum hmit can be taken (i) very easily by approaching the FM-FMD transition almost perpendicularly by tuning 
essentially one counterterm, and (ii) at a point on this transition line which is reasonably far away from the expansion 
point. This has been possible in ||7| and again in the present work. 

A central claim of the gauge fixing proposal is that it is universal, i.e., it should work with any lattice fermion action 
that has the correct classical continuum limit. This is because the central idea as discussed above is independent of the 
particular lattice fermion regularization. In the present paper we want to confirm the universality claim by applying 
the proposal to domain wall fermions with U{1) gauge group. For this purpose we have chosen the waveguide 
formulation pO| of the domain wall fermion and investigate in the reduced model. This model was investigated before 
without gauge fixing and the free domain wall spectrum was not obtained in the reduced limit |Tl| ]. Mirror chiral 
modes were found at the waveguide boundaries in addition to the chiral modes at the domain wall or anti-domain 
wall. 

In section II we present the gauge- fixed domain wall fermion action for a U{1) chiral gauge theory and then go to 
the so-called reduced model by switching off the transverse gauge coupling. In section III we perform a weak coupling 
pertu rbat ion th eory in the reduced model for the fermion propagators and mass matrix to 1-loop. However, in sections 
III C and [II E we have used special boundary conditions (instead of the actual Kaplan boundary conditions) to arrive 
at explicit expressions for the overlap of the opposite chiral modes. Our numerical results for the quenched phase 
diagram and chiral fermion propagators at the domain wall and anti-domain wall and at the waveguide boundaries 
are presented and compared with the perturbative results in section IV. We summarize in the concluding section V. 
In Appendix A, we describe the special boundary conditions used in sections HI C and [II D . In Appendix B we 
schematically discuss how using Kaplan boundary conditions one can arrive at the same qualitative conclusion about 
the 1-loop overlap of the opposite chiral modes. 



II. GAUGE-FIXED DOMAIN WALL ACTION 

Kaplan's free domain wall fermion action ^ on a 4 + 1-dimcnsional lattice is given by (lattice constant is taken to 
be unity throughout this paper), 

Sf = Y.^xW5-W5+M]^yi;Y (2.1) 

XY 

where ^ and are the fermion fields, and ^5 and W5 are respectively the 5-dimensional Dirac operator and the Wilson 
term, 

1 ^ 

i^5)xY = 2^^°" {Sx+&,Y - Sx-&,y) , 

Q = l 

5 

{w5)xy ^^Y1 (^x+&,y + Sx-&,Y - ^Sxy) , (2.2) 

Q = l 

The 7q's are the five hermitian euclidean gamma matrices, r is the Wilson parameter, X ~ {x, s), Y = (y, t) label the 
sites of the L'^Lg lattice and Lg is the extent of the 5th dimension: < s,t < — 1. We are interested in taking the 
continuum limit in the 4 space-time dimensions only. It is convenient to look at the 5-th dimension as a flavor space. 

With periodic boundary conditions in the 5th or s-direction {s,t = Ls^s,t — 0) and the domain wall mass M 
taken as 

M.XY = rn{s)SxY, where, (2.3) 



-mo, < s < Ls/2 
m{s) = 0, s = 0,L^/2 (2.4) 
Too, Ls/2 < s < Ls 

the model possesses a lefthanded (LH) chiral mode bound to the domain wall at s = and a righthanded (RH) chiral 
mode bound to the anti-domain wall at s = Ls/2. For m^Ls ^ 1, these modes have exponentially small overlap. 
The chiral modes exist for momenta p below a critical momentum pc, i.e. \fi\ < pc, where jP — 2 ~ cos(pp)] and 
= 4 — 2mo/r. Taking the Wilson parameter r = 1 the choice of toq is then restricted to < toq < 2. 
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A 4-diniensional gauge field which is same for all s-slices can be coupled to fermions only for a restricted number 
of s-slices around the anti-domain wall with a view to coupling only to the RH mode at the anti-domain wall. 
The gauge field is thus confined within a waveguide, 

WG = (s : So < s < si) 

with So = — 1, si = 1. (2.5) 

With this choice, {Ls — 2) has to be a multiple of 4. For convenience, the boundaries at (sq, so -I- 1) and (si, si -I- 1) 
are denoted waveguide boundary / and II respectively. 

The gauge transformations on the fermion fields are defined as follows: 

V'x ^ gx-iPl, i^l^ i'lgl s e WG 

r^^r,, HWG (2.6) 

where £ G, the gauge group. Other symmetries of the model remain the same as in [ pT| . 

Obviously, the hopping terms from sq to sq -I- 1 and that from si to si + 1 would break the local gauge invariance 
of the action. This is taken care of by gauge transforming the action and thereby picking up the pure gauge dof or 
a radially frozen scalar field (p (Stiickelberg field) at the waveguide boundary, leading to the gauge-invariant action 
(with (f,^ g^ipx and (|2^)) : 

Ip' {]p{U)-W{U)+m{s))i;' + ^ Ip' - w + m{s)) ^j' 



seWG s^WG 

E 



- y PlV"^^ +T^\PB.r") - y i^"' ^Pl^'^^ PRr') (2.7) 

where we have taken the Wilson parameter r = 1 and have suppressed all indices other than s. The projector Pl(r) 
is (1 =F 75)/2 and y is the Yukawa couphng introduced by hand at the waveguide boundaries. Ip (U) and W{U) are 
respectively the gauge covariant Dirac operator and the Wilson term in 4 space-time dimensions. ^ and w are the 
4-dimensional versions of ( |2.2|) . 

The gauge-fixed pure gauge action for U{1), where the ghosts are free and decoupled, is: 

Sb{U) - Sg{U) + SgfiU) + SctiU) (2.8) 

where, Sg is the usual Wilson plaquette action; the gauge fixing term Sgf (as proposed by Shamir and Golterman) 
and the gauge field mass counter term Set are given by (for a discussion of relevant counterterms see PJT^), 

\xyz X / 

S,,{U) = ~kY,{U^.x + uI), (2.10) 

Xfl 

where □(?/) is the covariant lattice laplacian and 

j^^^Y^fYtE^i+it^y (2.11) 



with = i (C/^, - C/t,) and ^ - 1/(2^^2). 

Sgf is not just a naive lattice transcription of the continuum covariant gauge fixing term, it has in addition 
appropriate irrelevant terms. As a result, Sgf has a unique absolute minimum at [/^j. = 1, validating weak coupling 
perturbation theory (WCPT) around g — or k = oo and in the naive continuum limit it reduces to ^ J d'^x{d^A^)'^ . 

Obviously, the action SsiU) is not gauge invariant. By giving it a gauge transformation the resulting action 
SBiy^lU^xipa:+ii) is gauge-invariant with U^^ ^ gxU^ 

x9x+fi and (fix 9x^x1 gx € By restricting to the trivial 

orbit, we arrive at the so-called reduced model action 
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Sreduced = Sf{U = 1) + Sb{^\ 1 Vx+jl) 

where Sf{U — 1) is obtained quite easily from eq.( |2.7| ) and 

X X 

now is a higher-derivative scalar field theory action. in ( ^.13| ) is same as in (2.11) with 



x+fj.t 



(2.12) 



(2.13) 



(2.14) 



In the following, we investigate the action (2.12) at y = 1 by analytical and numerical methods. Some numerical 
results with other values of y have been presented in |l3[ . The waveguide model strictly at y = would give rise to 
opposite chiral modes at the waveguide boundaries as can be seen from fermion current considerations Q (and also 
from numerical simulation) and would thereby spoil the chiral nature of the theory. It is an interesting question to 
investigate the model for < y < 1. Analysis of the results for small values of y{< 1) is tricky and will be discussed 
in a separate article [14| . 



III. WEAK COUPLING PERTURBATION THEORY IN THE REDUCED MODEL 

At y = 1, we carry out a WCPT in the coupling l/k for the fermion propagators to 1-loop. In order to develop 
perturbation theory, in reduced model, we expand. 



where, b = 1/V2k and 9^ is dimensionless, leading to 



S = Sf (V, ^; y) + 4"' (^) + 5^'"*^ (V', ^, 0\ y) 



where S'^^^'s are free actions and 5*^™*^ is the interaction part. 



A. Scalar propagator at tree level 



From Sg {6) one gets the free propagator for the compact scalar 6 , 

1 o K 



g{k) = 



F(fc2 +^2)- 



where, fc^ — 2sin(fcp/2). 



B. LL and RR fermion propagators at tree level and at 1-loop 



1. Tree level 



(3.1) 



(3.2) 



(3.3) 



With y Sp^-tjjj^p-jy — 1) is the free domain wall action (2.1). Free fermion propagators at y = 1 are obtained 



in momentum space for 4-spacetime dimensions while staying in the coordinate space for the 5th dimension following 
[ p5[ (results in |p5|-p7[ cannot be directly used because of difference in implementation of the domain wall (2.4)). The 
free action is written as, 



p.s,t 



(3.4) 



where, Mst = F{p)Ss,t + {Mo)st, (Afo)st = [l + m,{s)]6s,t-Ss+i,t, F{p) = I]^(l - cos(pp)), = sin(p^) and /J = 7^^^. 
The free fermion propagator can formally be written as. 



4 



A{p) = i^ + MPL+APPf 



= {-ii> + M^)PlGl{p) + {-ii> + M)PrGr{p) 



(3.5) 



where, 



Gl{p) 
Gr{p) 



E^pI + map 



Solution of Gl is obtained by writing ( p.6[ ) explicitly: 

+ l + B{sf] {GL)s.t - B{s + l){GL)s+i.t - B{s){Gl)s-u = 



(3.6) 
(3.7) 

(3.8) 



and similarly for Gr. In (3.8), B{s) — F{p) + l + m{s). We show only the calculations for obtaining Gl and henceforth 
drop the subscript L. 



Setting the notation as follows: 



G 

and G 



G- 
G^ 



B{s) = F{p) 
B{s) = F{p) 



mo 
mo 



a I 



for <s< Ls/2 - 1 
for Ls/2 < s <Ls~l, 



the equations for G are given by, 



if 



Git - a+G 



s+l,t 



a+G 



(3.9) 
(3.10) 



(3.11) 
(3.12) 



The ranges of s in eq s.(3.9, 3.1C ) for which G~ and G~^ are defined, ar e a pplicable only to the translationally 
invariant eqs.(3.11, 3.12). In gen eral for th e translationally noninvariant eq.( |3.8| ) we also define G" and G+ at s = 
and Ls/2, the ones excluded by ( 3.11 , 3.12 ). The ± superscript to G at s = 0, Ls/2 is decided by the translationally 
invariant s-sec tor fr om which s = or Ls/2 is approached in eq.( [3.§| ). We have used this notation for the boundary 
conditions eq.(3.15) below. 



The solutions of the eqs.(3.11, B.12) are expressed as sum of homogeneous and inhomogencous solutions: 



G± (p) = .9«(t)e-"±(^'> + g^i\t)e'^^(P> 



cosh[a±{p){\s-t\-l/2)] 
2a± sinh(Q!±(p)) sinh(a±(p)Z/2) ' 



where, I = Ls/2 and 



cosh(a±(p)) 



a± 



a± 



(3.13) 



(3.14) 



The third term in ( |3.13 ) is the inhomogencous solution. To avoid singularities in a±{p) when a± is zero further 
restricts the allowed range of toq to < toq < 1- In this paper we have taken toq = 0.5. 



In order to get the complete solution we need to determin e the unknown functions g')^^ (t) and g^^' [t) in (3.13), which 
are obtained by considering boundary conditions from eqs.( 3.11 , 3.12 ) at s = 0, 1, Ls/2 — 1, Ls/2, Ls/2 + 1, Lg — 1, 

aoGo,f(p) = a+Gl^^{p), 
(p2 + 1 + al)Go^tip) - a-Gltip) = So,t + aoG'l^_^ ^{p), 



,(2), 



it 



(3.15) 



with B[s) — F{p) + 1 = flo at s = 0, Ls/2. It is to be noted that these boundary conditions are significantly different 
from the ones given in ||l^ because of the difference in implementation of the domain wall. Gq ^ and Gj ^ (and the 
corresponding ones from Gr) are used to determine the free chiral propagators at the domain wall and anti-domain 
wall for comparison with numerical data in Figj^. We will see later that these chiral propagators do not receive any 
1-loop self-energy corrections. 
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S ubstit uting p 



eqs.( 3.15| ), we arrive at an equation of the form 



Fq and 2a± sinh(a±(p)) sinh(Q!±(p)//2) = X± and using the boundary conditions, 



A.g(i) =X(t), 



(3.16) 



where, g{t) — {g^_^^ g^'' g^^ 9+') is a 4-coniponent vector, X(t) is another 4-coniponent vector and A is a 4 x 4 



matrix as given below 
/ 



a-Ls/2 



Q+L3/2 



-floe 



-floe 



ao 

i^o-a_e-"- Fo-a_e"- -aoe-"+(^=-i) -ane"+ 



a+Ls/2 



and 



X(t) 



/ a+X+cosh[a+(|Ls - t| -1/2)] - aoX_ cosh[Q;_ (| ~ t\ -1/2)] 

aoX+cosh[a+(|L,/2-i| -1/2)] - a_X_ cosh[a_ (|L,/2 - t| - //2)] 

(5o.t + aoX+ cosh[a+(|Ls -l-t\- 1/2)] - FqX^ cosh[a_(| - t\ - l/2)] + 

a_X_cosh[a_(|l-i|-Z/2)] 

SLj2.t + a+X+ cosh[a+(|L,/2 + 1- t\ - 1/2)] - FoX+ cosh[a+(|L,/2 - t| - l/2)] + 

aoX- cosh[a_(|L,/2 - 1 - t\ - 1/2)] J 



\ 

The expHcit expressions for A and X(t) are obviously different from similar expressions given in [p"5|-pT| because of 
the differences in domain w all im plementation as already discussed earlier. 

The solution to the eqs.(3.16) is very complicated in general, particularly for finite Ls- However, g±(t) can be 
obtained for finite by solving the above equations numerically for different t values. This way we can easily 
construct the free fermion propagators at any given s-slice, including the zero mode propagators at s = i = 0, Ls/2. 

The solutions for (G r)s t and t he resul ting propagators are obtained i n ex actly the same way. However, in this case 
the explicit forms for (|3.8D, (3.13), (3.15) and matrices A and X(i) in (3.16) are obviously different. 



2. 1-loop 



Next we calculate the chiral fermion propagators to 1-loop. Half-circle diagrams which are diagonal in flavor space 
contributes to LL and RR propagator self-energies. However, the self-energies are nonzero only at the waveguide 
boundaries / and IL 

Retaining up to 0{b^) in the interaction term S*''"'^ {ip, ^,d]y = 1) in (3.2), we find the vertices necessary to calculate 
the self-energies to 1-loop. 



' R 



R > 



p, So + 1 p- k, So p. So -I- 1 
FIG. 1. 1-loop self-energy contribution to LL propagator at WG boundary /. 

The LL propagator on the (so + l)-th slice at the waveguide boundary / receives a nonzero self-energy contribution 
from the half-circle diagram. 



SL(P)),, = (04^' [-lj,{p- k),PLGL{p- k)]^^^^^^g{k) 6s,so + lS, 



BZ (27r) 

^,T.[^^kp-k) 



't,So + l 



So, So fc2(^^2 _|_ ^2^ 



s,so + l'^t,so-l-l 



(3.17) 
(3.18) 
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where the expression in the square bracket in ( 3.17 ) is the free RR propagator [iSj^pJap^so on the so-shce. Eq.(3.17 



assumes i nfini te 4 space-time volume while in eq.( |3.18 ) a finite space-time volume L is considered 



Using ( |3.18 ) we numerically evaluate the analytic 1-loop propagator on a given finite lattice, using 



'Li 



= 5 



(0) 
LL 



S 



(0) 
LL 



^LL 



s 



(0) 
LL' 



(3.19) 



and in the section IV compare with nonperturbative numerical results. To avoid the infra-red proble m in t he scalar 
propagator, we use anti-periodic boundary condition in one of the space-time directions in evaluating ( 3.18 ). 

In a similar way, 1-loop corrected RR or LL propagators are obtained at all the s-slices of the waveguide boundaries 
/ and //, i.e., at the slices so, so + 1, si and si -I- 1. 



C. Fermion mass matrix at tree level and 1-loop 

Another issue of interest is the spread of the wavefunctions of the two chiral zero mode solutions along the discrete 
s-direction and their possible overlap. A finite overlap would mean an induced Dirac mass. The extra dimension, as 
already pointed out in the discussion following ( ^.2|) , can be interpreted as a flavor space with one LH chiral fermion, 
one RH chiral fermion and (Ls/2 — 1) heavy fermions on each sector ofl<s<Ls/2 — 1 and Ls/2 + l<s<Ls — 1- 



1. Tree level 



For the spread of the zero modes at the tree level, one needs only to solve MqUj^ = and M^un = [plp^] 
where Mq — M{p — 0) (keeping the momentum p non-zero is unnecessary in this discussion). However, for radiative 
correction on the domain wall mass m{s), the heavy mode spreads are also needed. Accordingly we consider flavor 
diagonalization of MqMq {Mq is not hermitian) as in [ p6| , p7[ : 





[Af] 




).,(*"'), ^ 







(3.20) 
(3.21) 

The index j for the eigenvalues and the eigenvectors is basically a flavor index, but unlike s and t which vary from 
to is — 1, it is taken symmetric around the domain wall j = 0. Index j varies from —Ls/2 to Ls/2. From periodic 
boundary condition, Ls/2 and —Ls/2 are the same point in flavor space corresponding to the anti-domain wall. It is 
to be noted that j appears explicitly in the heavy mode solutions below. However, it may be pointed out that the 
index j need not be chosen this way. One could also define j in the same way as the flavor indices s and t, only in 
that case the explicit solutions below would have a different appearance and to our taste less tractable. 

To solve the above tree level eigenequations (later also at 1-loop), we should ideally use the Kaplan boundary 
conditions because we used in our action the Kaplan way of implementing the domain wall. However, that would 
make an explicit calculation of the heavy eigenmodes quite complicated and almost intractable. We have hence devised 
a special set of boundary conditions (look at Appendix A for details on the boundary conditions) with the property 
that no information can be passed through the wall and the antiwall. This makes the calculation less cumbersome and 
explicit expressions can be obtained in manageable forms. We stress that these are not the actual boundary conditions 
in the particular domain wall implementation we have taken. However, we show in Appendix B that using the correct 
boundary conditions {i.e., the Kaplan boundary conditions) also would lead to the same qualitative conclusions about 
the nature of 1-loop corrections to the eigenvalues, stability of the zero modes, and particularly the 1-loop overlap of 
the opposite chiral zero modes. 

With our special boundary conditions (used only in sections [II C and HID), the explicit form of the eigenfunctions 
are obviously different from (l^Jl^ and are obtained as. 



Zero (L-handed) mode: 

(j = 0) 
(j = Ls/2) 

with. 



= A exp(-as) < s < Ls/2 

S 

= A Gxp [—a{Ls — s)] Ls/2 < s < Ls 

/ 1 - exp(-2a) 
\1 - exp{~aLs) 
exp(— a) = 1 — mo 



, 1/2 



(3.22) 
(3.23) 

(3.24) 
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Heavy mode ^ 0, Ls/2): 

(0<j<L,/2) (^f^) = y^sin/3j(s-l) < s < i,/2 (3.25) 

(-L,/2<j<0) = W— sin/3_j(L,/2-s + l) Ls/2<s<L,. (3.26) 

V -tJ.s 



The parameters involved are as follows: 



[xf^^'^ = 1 + a(s)2 - 2a(s) cos (3.27) 
a{s) = l + m{s) (3.28) 



and for nontriv ial he avy mode solutions, (3j — 2T:j / with \j\ ^ 0, Ls/2. 
Solutions to ( 3.21 ) for all j and s easily follow: 



It is obvious that the eigenvectors <t>flii and 'I'j'^o correspond respectively to the LH and RH chiral zero modes at 
the domain wall and the anti-domain wall in the region < s < Ls/2. Similarly the eigenvectors <P^^1]^ y2 ^'^'^ *^j=i /2 
correspond to the same chiral zero modes at the domain wall and the anti-domain wall in the region Ls/2 < s < Lg. 
The |j| 7^ 0, Ls/2 eigenvectors arc for the heavy flavor modes. All solutions are real. 

The overlap of the opposite chiral modes in the region < s < Ls/2 does not depend explicitly on s and is given 
at the tree level by, 

1>fJo -^fio = exp(-ciLs/2) (3.30) 

and similarly in the region Ls/2 < s < L^. The exponentially damped mixing of the LH and RH chiral modes does 
not induce any Dirac mass at the tree level for large Lg. 

It is also noted that (p^^^ and C'^"^ diagonalize Mq and A/q in the following manner: 

(^°0s (^?). - ^r^..' - K),, (^?). (^-^i) 



The way the formalism is set up, the absolute sign of A^*^^ in eq.(3.31) is arbitrary. 



2. 1-loop 

Flavor off-diagonal tadpole diagrams produce the self-energies for the LR and RL parts of the fermion propagator. 
Again, the self-energies are nonzero only at the waveguide boundaries / and //. 



L ^ \_ _/ ^ R 

p. So + 1 So 
FIG. 2. Tadpole contribution to 1-loop LR propagator at WG boundary 7. 

For the LR propagator connecting sq and so -|- 1 at the waveguide boundary /, the self-energy contribution from 
the tadpole diagram is given by, 

= ^b^PLTSs^g,St,so+i (3.33) 
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where T ^ 0.04 is the tadpole loop integral. 

Similarly the self-energy contribution to the LR propagator at the waveguide boundary // connecting si and si + 1 
comes from a tadpole diagram and is given by, 



The mass parameter Mq gets modified at 1-loop as: 



(3.34) 



{Mo)stPL {Mn)stPL = [Mo)stPL + [-(SiflXO)).*] + [-(Sf«(0)),t] 
= {M^)stPL + h\i:LR)stPL- 



(3.35) 



^LL ~ ^ identically. {Ml)stPR gets modified accordingly. 

Because of our use of the special boundary conditions (Appendix A) which do not allow any communication 
through the walls, 1-loop correction to the eigenvectors and eigenvalues involves either sq (waveguide boundary /) or 
si (waveguide boundary //) depending on which sector of j is considered. In Appendix B we show that if the actual 
(Kaplan) boundary conditions are used, one gets contribution at 1-loop level also from another diagram, called there 
the global loop diagram. However, the conclusions are qualitatively the same. 



D. Mass matrix diagonalization at 1-loop 



At 1-loop level we organize the corrections to the eigenvectors and the eigenvalues of the fermion mass matrix 
squared as follows: 



(/.(°)^0=(l + 6Vi))(/.(°), 
= (l + 62$(i))$(o), 



A^^^ is found to be 



A« = r*?) 



The eq.(3.31) gets modified in 1-loop as 

(<i>,)s(Mo)st (</>/)* = \,5,,r+0{h^) = (</>,). (M]).t(<i>,Ot. 



The 1-loop correction to the eigenvalues, {5\)j — , is given by: 

{6\), = — r sin/3j(L,/2 - - 1) sin/3,so + 0{b^) 

952 

= rsin/3j(si -I- 1) sin/3j (L,/2 - si) + 0(6*) 



yt2rcxp[-5(L,/2 + l)] +0(6-*) 



Ls — 'oo 







for l<j< (is/2-1), 
for -(L,/2-l)<j<-l, 
for J = 0, L,/2. 



(3.36) 
(3.37) 
(3.38) 

(3.39) 
(3.40) 
(3.41) 



We notice in eq.(3.41) that the 1-loop correction to zero mode eigenvalue is exponentially damped and the zero modes 
are hence perturbatively stable. 

The 1-loop expression for the chiral zero mode at the domain wall, in the region < s < Ls/2, is. 



t'j=o)s = A 



is/2-1 



2b^ 1 

cxp(-(5s) - — Texp(-a(so + 1)) ^ -7^ sin/3jv(Lj,/2 - so - 1) sin/3j/(s - 1) 



^ x(o) 



(3.42) 



We obtain a similar 1-loop expression for the chiral mode at the anti-domain wall. 

The overlap of the opposite chiral modes in the region < s < Ls/2 at 1-loop is given by. 
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<i>j=o 0j=o = A'^ exp(-<5L,/2) - ( — ) A^T Ti cxp[-a(so + 1)] 
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A^TT2 cxp[-5(L,/2 - so)], (3.43) 



where, 



= '^-^e-xp[-a{Ls/2 - s)] sin /3jv(s - l)sv[i(3.y{Ls/2 - sq - 1) 



= ^ — exp(— as) sin/3j'So sin/?^' (is/2 — s — 1). 



Eq.(3.43) clearly shows that 1-loop corrections to the mixing of the LH and RH modes are also exponentially damped. 
This guards against any induced Dirac mass in the domain wall model for large enough Ls and the waveguide 
boundaries / and // chosen approximately equidistant from the domain wall and the anti-domain wall. In this 
context we point out that in the Smit-Swift model a shift symmetry for the singlet fermion ensures that no fermion 
mass counter term is needed. 



Using the actual Kaplan boundary conditions would make the relatively nice explicit expressions in eqs.( ^j.3£ - 3.43| ) 



much less manageable, to say the least. Diagonalization of the fermion mass matrix upto 1-loop using the Kaplan 
boundary condition is discussed in Appendix B. 



IV. NUMERICAL RESULTS 

In the quenched approximation, we have first numerically confirmed the phase diagram in ||l9|| of the reduced model 
in (k, k) plane. The phase diagram shown schematically in Figj^ has the interesting feature that for large enough k, 
there is a continuous phase transition between the broken phases FM and FMD. FMD phase is characterized by loss of 
rotational invariance and the continuum limit is to be taken from the FM side of the transition. In the full theory with 
gauge fields, the gauge symmetry reappears at this transition and the gauge boson mass vanishes, but the longitudinal 
gauge dof remain decoupled. In Fig.^ PM is the symmetric phase and AM is the broken anti-ferromagnetic phase. 
The numerical details involved in reconstruction of the phase diagram and the fermionic measurements that follow 
will be available in [Q. 



\ 


FM 


• 




\ PM \ 


FMD 


_ -2 






_ -4 








AM \ 






\ 0.1 


0.2 0.3 



K 

FIG. 3. Schematic quenched phase diagram. 



For calculating the fermion propagators, as in |7j we have chosen the point n — 0.05, k — 0.2 (gray blob in FigJ^). 
Although this point is far away from k = oo, around which we did our perturbation theory in the previous section, 
the important issue here is to choose a point near the FM-FMD transition and away from the FM-PM transition. The 
results below show that for the fermion propagators there is excellent agreement between numerical results obtained 
at K = 0.05, k = 0.2 and perturbation theory. 

Numerically on 4^16 and 6'^16 lattices with Lg — 22 and mg — 0.5 we look for chiral modes at the domain wall 
(s = 0), the anti-domain wall (s = 11), and at the waveguide boundaries (s = 5, 6 and s = 16, 17). Error bars in all 
the figures are smaller than the symbols. 
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Fig.^ shows the RR propagator ISrrI and the LL propagator \Sll\ at the domain and anti-domain wall as a 
function of a component of momentum p4 for both p = (0, 0, 0) (physical mode) and (0, 0, tt) (first doubler mode) at 
y = 1. From the figures, it is clear that the doubler does not exist, only the physical RR {LL) propagator seems to 
have a pole at p ~ (0, 0, 0, 0) at the anti-domain (domain) wall. In all the figures, NS, FT and FF respectively indicate 
data from numerical simulation, from perturbation theory and from free fermion propagator by direct inversion of the 
free fermion matrix. 



3.0 
2.5 

2.0 

\Srr\ 
1.5 

\Sll\ 
1.0 

0.5 





A 



\Sll\ 



A 
□ 



domain wall (s = 0) 

6^ X 16, (p = 0,0,0) FF 
8^ X 20, {p = 0,0,0) FT 
6^ X 16, (p = 0,0,0) NS 
6^ X 16, (p = 0,0, tt) NS 

(k, k) = (0.05, 0.2) 



E] 



\Srr\ 



□ . 



0.4 



0.8 



1.2 



• , ■ * 
JS 2 



P-l 



3.0 
2.5 

2.0 

\Srr\ 
1.5 

\Sll\ 
1.0 

0.5 





\Srr\ 



antidomain wall (s = 11) 

A 63 X 16, (p = 0,0,0) FF 
□ 8^ X 20, (p = 0,0,0) FT 

• 6^ X 16, (p = 0,0,0) NS 

* 6^ X 16, (p = 0,0, tt) NS 

(k, k) = (0.05, 0.2) 



E] 



\Sll\ 



□ . 



A 



■ ★ 



0.4 0.8 1.2 1. 



2.0 



FIG. 4. Chiral propagators at domain wall 5 = and at anti-domain wall s = 11 {La = 22; a.p.b.c. in L4, y = 1.0). 
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FIG. 5. RR propagator at waveguide boundary s 
in Li, y = 1.0). 



5 and LL propagator at waveguide boundary s = 6 {Ls = 22; a.p.b.c. 



For Fig.^, PT also means zerot h order perturbation theory, i.e., numerical solution of propagator following eg. ( 3. If ) 
(as noted before in subsubsection III B 2 the self-energy contributions to the LL and RR propagators are nonzero only 
at the waveguide boundaries, the propagators in Fig.^ do not get any 1-loop correction). We have PT results also 
for 6^ X 16 lattice but have chosen not to show them because they fall right on top of the numerical data. Instead 
PT results are shown for 8'^ x 20 lattice for which the pi points are distinct. The dotted line in all figures refer to 
the propagator from PT using a 256'^ x 1024 lattice. The curves stay the same irrespective of methods or lattice size. 
Based on the above, we can conclude that there are only free RH fermions at the anti-domain wall, and at the domain 
wall there are only free LH fermions. 

Fig.p show no evidence of a chiral mode at the waveguide boundaries s = 5 and 6 and excellent agreement with 1-loop 
perturbation theory. Here too doublers do not exist. Actually the agreement with the FF method (direct inversion 
of the free domain wall fermion matrix on a given finite lattice) is also excellent, because the 1-loop corrections are 
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almost insignificant. For clarity, in Fig.|5| we have not shown the LL propagator on s = 5 and the RR propagator at 
s = 6, but conclusions are the same. 

Similar investigation at the other waveguide boundary s = 16, 17 also does not show any chiral modes. Previous 
investigations of the domain wall waveguide model without gauge fixing [ pl have shown that the waveguide boundaries 
are the most likely places to have the unwanted mirror modes. This is why we have mostly concentrated in showing 
that there are no mirror chiral modes at these boundaries, although we have looked for chiral modes everywhere along 
the flavor dimension. In fact, we do not see any evidence of a chiral mode anywhere other than at the domain wall 
and the anti-domain wall. 



V. DISCUSSION 

We have followed the gauge fixing proposal of Shamir and Golterman and applied it to domain wall fermions for a 
U{1) LxGT. By switching off the transverse gauge dof, we arrive at the so-called reduced model. We have determined 
the quenched phase diagram of the model and confirmed that there is a continuous phase transition from a broken 
symmetry ferromagnetic (FM) phase to a broken symmetry rotationally noninvariant ferromagnetic directional (FMD) 
phase with the properties that at this transition the longitudinal gauge dof, i.e., the if fields get decoupled. We have 
come to this conclusion by performing a WCPT for the fermion propagators and comparing them to nonperturbative 
numerical simulations. 

Let us now contrast this with the previous attempt jTll of lattice regularizing a chiral gauge theory with domain 
wall fermions. No gauge fixing was done and in the reduced model, a combination of analytic and numerical methods 
showed that mirror chiral modes were dynamically generated in the theory. 

As in the Smit-Swift model ||l|, the mirror chiral modes are a reflection of the undesired presence of the longitudinal 
gauge dof in the continuum limit. Without gauge fixing, these radially frozen group- valued scalar ip fields are generally 
nonperturbative or rough for any value of the transverse gauge coupling, even in the reduced model limit. If the 
continuum limit is taken at the FM to a symmetric paramagnetic (FM) phase transition, the if fields survive with 
radial modes and physical effects of them coupling with the rest of the theory become manifest in the form of mirror 
chiral modes etc.. 

In the gauge fixing approach of Shamir and Golterman, care has been exercised so that the gauge fixing term is 
not just a naive lattice transcription of the continuum gauge fixing condition. There are appropriate additions of 
irrelevant terms in the covariant gauge fixing term so that a unique perturbative vacuum exists. With this gauge 
fixing in the reduced model limit the ip fields become smooth and can be perturbatively expanded as 1 + 0{1/\/k) + ■ ■ ■ 
around k = oo. We have found in our investigation that as long as k is taken sufficiently large so that a continuum 
limit exists from the FM phase to the FMD phase (not the PM phase), the model is as good as at the perturbative 
limit K = oo. This seems to hold true in our particular implementation, i.e., the domain wall fermion case, in a strong 
sense, because the perturbative 1-loop corrections (i.e., the leading nontrivial corrections) to the fermion propagators 
are found to be negligible (Figs. 4 and 5). In gauge fixing the Smit-Swift model, however, 1-loop corrections to the 
fermion propagators were small but not negligible 0. In the present investigation with domain wall fermions, to the 
accuracy of our calculations, the reduced model spectrum is that of a free domain wall model with absolutely no trace 
of the (f fields. 

In our investigation of the reduced model, only a counterterm with coefficient k was sufficient to reach the FM-FMD 
phase transition by decreasing k from the FM-side. A dimension-three fermion mass counterterm was not needed, 
because domain wall fermions have the robust property that mixing between the opposite chiral modes at the domain 
wall and the anti-domain wall is exponentially damped and is negligible for large L^. 

We have carried out and presented the perturbative results for fermion propagators and mass matrix in reasonable 
detail, because, although the technique employed is not new, the explicit results are available mostly for only the 
QCD (wall) implementation rather than the wall-antiwall implementation suitable for a chiral gauge theory. With 
transverse gauge fields back on (and fermions in an anomaly-free representation), the full gauge-fixed domain wall 
model of LxGT is ready for a perturbative treatment with the same techniques as used in this paper if the transverse 
gauge coupling is perturbative. The model is also ready for a nonperturbative treatment by numerical simulation for 
strong transverse gauge coupling. 

Our numerical computations have been done in the quenched approximation and they agree for the fermion prop- 
agators very well with zeroth order perturbation theory with the 1-loop corrections almost negligible in our case. 
Effects of fermion loops would enter these calculations at least at the 2-loop level. Inclusion of dynamical fermions in 
our nonperturbative numerical investigation should not affect our results about the spectrum of the model at all if 
the relevant part of the phase diagram remains qualitatively the same. 
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So far the gauge fixing method is applicable only to the abelian theory. Extension to a nonabelian gauge group is 
nontrivial and is being pursued at the current time § . 

As commented at the end of section II, the reduced gauge fixed domain wall model with < y < 1 is interesting 
because at y — the model is known to have mirror modes at the waveguide boundaries. This will be taken up in a 
separate publication jl^ . 
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APPENDIX A: 



In this appendix, we describe the boundary conditions used to flavor-diagonalize eqs.( 0.2C| ) and (3.21) to obtain 



the results presented in subsection III C and in III D , The boundary conditions make sure that no information passes 
through the domain wall and the anti-domain wall. This means that the eigensolutions corresponding to a given 
source flavor index s belonging to a particular segment < s < Ls/2 ot Ls/2 < s < Lg will be restricted to 
that particular segment. As already stated before, these are not the actual boundary conditions under which the 
numerical simulations or the analytic calculation for the LL/RR chiral propagators have been performed, for the 
flavor diagonalization this enables us to calculate the heavy mode solutions and overlap of the chiral modes explicitly. 
In appendix B we have schematically described how the calculations would go in case of the actual Kaplan boundary 
conditions, but to obtai n explicit expressions would be a horrendous task. 

We consider eg. ( 3. 20 ) in the region < s,t < Ls/2. Obviously the index j is also in the same region. Writing 
eq.(3.2C) at s = Ls/2 gives (by force using the general expression for MqAIq away from the boundaries of the region). 



(Al) 



Again writing eq.( t3.20 ) at s = Ls/2 gives (this time using the expression for MqMq at the boundary of the region 
and dropping terms that take the solutions across the anti-domain wall), 



+ a,L./2,^-(Af)]«')^^^^-a,.,/2-i)(*f)^^^^__ 

For both the eqs.(Al) and ( [A^ ) to be true, the following must be true: 



= 0. 



(A2) 



(A3) 



which is a boundary condition. 



Two equations corresponding to the eqs.(Al, A2), this time at s = 0, are: 



(A4) 



Hence we get the boundary condition at s = to be 



= 0. 



(A5) 



(A6) 



For the same eigenequation eq.(3.20) in the other region Ls/2 < s,t < Ls, we similarly get two more boundary 
conditions as follows: 
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«')^^^^-a(W2-i)(*f),,„_,=». 



= 0. 



(A7) 
(A8) 



The se b ounda ry c onditions are similar to those given by Shamir pU, except that we get two sets of them, i.e. (|A3|), 



( A6), (A7) and (A8), cor respo nding to two sectors of s. Another couple of sets of boundary conditions are obtained 
for the eigenequation eq.(3.21) in the two regions in the same way. 



APPENDIX B: 



Here we schematically show that using the correct boundary conditions (in this case the Kaplan boundary conditions, 
instead of the special boundary conditions explained in Appendix A) does not qualitatively change our conclusions 
about fermion mass counterterms. 



The use of special boundary boundary conditions in sec. Ill C meant that the 5-th or flavor dimension was basically 
divided by the domain and the antidomain wall into two segments which did not communicate to each other. As a 
result a wave incident on a wall got fully reflected. In general, there will be transmission to the other side of the wall. 

The chiral zero modes are the same as before (except for a tiny change in the normalization) except that they are 
now allowed to decay all the way around the s-space. 

Let us assume then the following trial eigenfunctions (f>^ and $^ for the heavy modes: 



-±m 



A±j exp(— i/3±js) + B±j exp(j/3±js) 



($^- )s = C±j exp(-j/3±js) + D±j exp(i/3±js) 



(Bl) 
(B2) 



The subscripts ± go with the regions in s-space with corresponding signs of m(s) in eq.(2.4). The wavevectors f3 are 
obtained by trying these trial solutions in the eigenvalue eqs. ( 3.20| ), ( |3.21 ) and are given as: 



2a(s) cos /3±j = 1 + d{sf - 



(B3) 



We proceed essentially the same way as in the case with the free propagators except for the fact that here we do 
not have the inhomogeneous part of the solutions. The boundary conditions used to determine A±j, B±j, C±j, D±j 
are (omitting the subscript j and also the tree level indicating superscript (0) for convenience): 



(1 



(1 



;;2 



ao(,' 



aocp 



= 



'LJ2-1 



(B4) 
(B5) 
(B6) 
(B7) 



(1 



(1 + ao - X^)% - ao$J,_i - a-$r = 



A^X/2 



is/2-1 



«+*t,/2+l 



(B8) 
(B9) 
(BIO) 
(Bll) 



where a± = 1 ± mo and ao = = 1 from eq.( 3.28 ). 

This leads to two sets of linear homogeneous equations for the amplitudes: 



PA 

Q-e 



where P and Q are 4x4 matrices and A ~ (^-, A+, B+) and 8 — (C_ 

solutions, det P = and det Q = leading to (reviving the subscript j) 



(B12) 

D_, C+, -D+). For nontrivial 



R , 
P±3 = 1-3 



(B13) 



14 



Solutions to the eqs.(B12) are of the following general form (dropping the subscripts j again): 



B_=/Pa_, A+ = /fA_, i3+ = /fA_ (B14) 
I?_=/1j'^C_, C+-/fC7-, D+ = /^^C- (B15) 

where fj^'^'^"^ are complex numbers with finite magnitude (order 1). and C_ are determined from normalizations 
of the respective eigenfunctions and would lead to a 0{1/Ls) factor. 

At 1-loop with the Kaplan boundary conditions, there is also a contribution to the fermion self-energy for the LR 
and RL parts coming from a flavor off-diagonal half-circle diagram (we shall call it a global-loop diagram) where the 
scalar field goes around the flavor space connecting fermions at the waveguide boundaries I and II. The global-loop 
diagram originates from the fact that the ip field that couples the fermions at the waveguide boundary / is the same 
if field coupling the fermions at the waveguide boundary //. 



• R L > R 



p, so + 1 So p-k si + 1 P, si 

FIG. 6. Global loop contribution to 1-loop LR propagator connecting WG boundaries / and //. 

Self-energy contribution from the global-loop diagram for the LR propagator is: 

- (Si'fl(p))^^ - b'PL J^^ ^ [Aft(p_ fc)Gi(p- fc)],^+,_^^ g{k) Ss,.Jt,so+i (B16) 
^b^PLn5s,,,St.,s,+i (B17) 



where TZ is the loop integral in eq.( B16D 



Now the calculations for the 1-loop corrected eigenvalue proceeds exactly as in subsection [IIP , except that it now 
includes also the contribution from the global-loop diagram: 



where, for \j\ ^ 0, Ls/2, 



A,=Af (Si«)^.^. (B18) 



{•SLB:)n^{<ff')s{^LB)st{cf>f'), (B19) 



(($7<"').„ (07'"').„,, + (a>;<"')., (0;<"')..+r) +mr)s. (07")^o.i, (B21) 



and similarly for the chiral zero modes. 

Using the expressions for the eigenfunctions at the specific values sq, si, sq + 1 and si + 1, wc arrive at the 1-loop 
correction to the eigenvalues: 

(a), ^ I (rji«(.oJi^/i^^) + rfif (.lJi^/F^/i,^^/r) + 7^^f , \j\^o, ^ (B22) 

^b'^Texp[-dt{Ls/2+l)]+b'^nexp[~dt{so + l)-d{Ls/2~si)], j = 0, y (B23) 

The details of the functions n = 1,2,3 are not illuminating. Correction to the zero mode eigenvalue clearly 

shows the exponential damping and for large Ls it is negligible. 
The 1-loop chiral zero mode wavefunction at the domain wall is, 
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~ exp(— as) 



A 



(0) 



(0)^ 



exp(— as) 

__r exp[-ci(so + 1)] ei(s) - expha(si + 1)] U^) ~ ^>'7^expha(so + 1)] ^(s), 



(B24) 
(B25) 

(B26) 



where the weU behaved functions ^i, i — 1,2,3 contain the information of the heavy mode wavefunctions and need not 
be written down expHcitly (it is to be noted j ^ in eq.(B25)). The 1-loop correction to the zero- mode wavefunction 
is also explicitly damped and for large negligible. 

The 1-loop correction to zero-mode wavefunction at the anti domain wall is similarly found to be damped. As a 
result the overlap of the opposite chiral zero modes at the 1-loop level is also exponentially damped. 
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